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We present the results of an extended scaling test of quenched Wilson twisted mass QCD. We 
fix the twist angle by using two definitions of the critical mass, the first obtained by requiring 
the vanishing of the pseudoscalar meson mass tops for standard Wilson fermions and the sec- 
ond by requiring restoration of parity at non-zero value of the twisted mass ji and subsequently 
extrapolating to ji — ► 0. Depending on the choice of the critical mass we simulate at values of 
P G [5.7,6.45], for a range of pseudoscalar meson masses 250 MeV < mps < 1 GeV and we 
perform the continuum limit for the pseudoscalar meson decay constant /ps and various hadron 
masses (vector meson my, baryon octet m oct and baryon decuplet m^c) at fixed value of romps. 
For both definitions of the critical mass, lattice artifacts are consistent with O(a) improvement. 
However, with the second definition, large 0(a 2 ) discretization errors present at small quark mass 
with the first definition are strongly suppressed. The results in the continuum limit are in very 
good agreement with those from the Alpha and CP-PACS Collaborations. 
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1. Introduction 

Twisted mass QCD (tmQCD), whose action reads 

S[U,Y,¥] =a A Y,W(x)(Dw+m + inY 5 T 3 )Y(x) , (1.1) 

X 

has been proposed as an alternative to Wilson QCD because it is not affected by the problem 
of unphysical zero modes and can lead to simplifications of the operator mixing pattern [jl]|. At 
maximal twist (i.e. by setting mo to its critical value m c up to 0(a), with p now the bare quark 
mass) it has been shown [^] that parity even correlators (and thus energies and matrix elements) are 
automatically 0(a) improved. Due to these properties, tmQCD is a very interesting candidate for 
dynamical simulations at small quark masses [^Jj. However, right at small quark masses p, <a , 
large lattice artifacts have been observed when using a definition of m c obtained by requiring the 
vanishing of mp S with standard Wilson fermions (in the following we will call it m^ mn ) [f|]. This 
kind of lattice artifacts are obviously affecting also dynamical simulations. An analysis based on 
Wilson Chiral Perturbation Theory (W^PT) [Q] suggested a definition of m c suitable to reach quark 
masses p ~ a 2 . This definition is obtained by requiring the vanishing of the PC AC quark mass 

mpCAC= ^ 2I x (^(^(0)) (L2) 
where P a = \fry 5 ^y,A^ = \j/y^ y 5 f y. 

An analysis a la Symanzik beyond 0(a) shows the presence of 0(a 2 ) cutoff effects enhanced 
at small pion mass, the most dangerous of which are of order (a/rrip S ) 2k k > 1 which signal the 
presence of possibly large lattice artifacts in the mjs s < a regime [Q] (see also The result of 
this analysis is that there are two ways of reducing these large lattice artifacts: fixing the 0(a) 
ambiguity in m c as proposed in Refs. ^ (we will call this definition m^ CAC ) or add the clover 
term with non-perturbatively determined csw coefficient (and the corresponding value of m c ). We 
investigate here the first proposal (see Refs. [||] for the second possibility). 

2. Determination of m^ CAC 

'Mpcac depends, in the neighborhood of a given estimate of m c (e.g. m? lon ), smoothly on 
both mo and p. Moreover, in the quenched approximation, a multiple mass solver can be used to 
compute the fermion propagator for different p values at a given value of mo- In view of these two 
facts, the procedure we adopted to determine m^ CAC is the following [gj: 

1. choose n values of pL (with p. > a 2 ) and n' values of mo (in the vicinity of m V c° n ) which cover 
the range where mpcAc(p) is close to zero (in the present case n = 9 and n' = 4). 

2. find the values of m c (p) at which mpcAc(M) i s zero ( see Fig- |} a ) 

3. extrapolate m c (p) to p = (see Fig. |[a) 2 . 

In Ref. [|6|] one can find a theoretical analysis that justifies this procedure and shows that the value 
of m^. CAC in which we are interested can be obtained by a linear extrapolation from the region 
p > a 1 down to p = 0. m^ CAC has been determined for j3 G {5.7,5.85,6.0,6.2} using statistics of 
O(100)-O(200) gauge configurations. 

'powers of Aqcd required to match physical dimensions are in the following understood 
2 in Ref. po| ] for each simulated value of the corresponding value of m c (jli) were used. 
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Figure 1: a. Determination of m c (ji) {k c = 2am c + 8) for given values of /i and j3 and extrapolation of 
m c (n) to ji = 0. b. fps as function of nip S at j3 = 6.0 for the two definitions of m c (mf on and ;«^ CAC ). 



3. Chiral behaviour at fixed a 

The effectiveness of in reducing the large lattice artifacts observed at small quark mass 

is evident when considering the chiral behaviour of two simple observables: fp$ and mps. fps can 
be extracted, without need of renormalization constants, by using the exact lattice PCVC relation 
( d ^(x)O(0)) = -2ne 3ab (P b {x)O{0)} (where d* is the lattice backward derivative and V™ the 
point-splitted vector current). It follows that, at maximal twist, 

M = ^\(0\P a \PS)\. (3.1) 

In Fig. |].b one can compare the chiral behaviour of fps obtained with either mf on or m^ CAC . One 
sees immediately that m^ CAC reduces the large lattice artifacts present at small /J. when using m^ lon 
(fps is predicted to be linear in mp s at one loop in quenched #PT). 

Using the integrated PCVC relation, it is also possible to prove ^ that, by using m^ CAC and 
in the region jx > a 2 , m 2 PS is linear with ji. up to small a 4 cut-off effects 3 (i.e. 0(a 2 ) lattice artifacts 
are proportional to ji). This is qualitatively confirmed by our data as shown in Fig. |2|a. 

4. Scaling behaviour 

We present now results concerning the scaling behaviour of /ps?"o, wxyro, m oct ro and m^cH) at 
fixed value of mpsrrj (where ro is the Sommer scale). More details concerning meson quantities 
can be found in Ref. [jTl||. my has been extracted by using either the spatial component of the axial 
vector or the temporal component of the tensor as interpolating operators (in the following we will 
quote only the latter, which systematically present lower statistical errors). m oct and have been 
extracted by using respectively e ABC ((d A ) T Cj5U B )u < ^ and e ABC ((u A ) T CYkU B )u^, as interpolating 
operators. The parameters of the simulations can be found in Tab. ||. We have simulated quark 
masses pL corresponding to 235MeV < mps < l.OGeV (where the scale ro, as will be explained 
below, has been fixed through the p mass). The scaling behaviour of fpsro, as shown in Fig. [|b, is 
clearly linear in (a/ro) 2 . However, m V c° a gives large 0(a 2 ) effects at small masses, effects that are 

3 chiral logs and other 0(rap S ) contributions are here assumed to be negligible 
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drastically reduced by the use of rriF c c . This obviously influences the scaling region which starts 
at j3 = 6.0 for m^ lon and at j8 = 5.85 for m^. CAC . We perform indipendent continuum extrapolations 
for the two choices of m c and the results are in good agreement. In the case of «i? lon , due to the 
highest slope for the lowest quark mass, we needed an additional point at j8 = 6.45 in order to 
control the extrapolation, myro, m oc tro and m<i ec ro have been thus computed by using only m% . 

Since we have simulated down to m PS of 235 MeV, finite size effects (FSE) can be quite 
relevant. In order to check for FSE we performed two additional simulations at j8 = 5.85 on 



volumes of 12 3 x 32 and 14 3 x 32 in order to extend the results of Ref. [13] at smaller masses. For 
meson quantities (mps, fps and my) FSE are negligible for all the quark masses starting from the 
third smallest one; on the two smallest masses they are in practice below the statistical accuracy 
of our data. For baryon masses, instead, FSE are very large for the two smallest masses and still 
relevant on the next three. Since the sensitivity required to study FSE on the smallest two masses 
is computationally very expensive, we chose here to correct only the data from the third value of jl 
on (corresponding to nips = 375 MeV). The results for the scaling behaviour of m oct ro are shown 
in Fig. [|.a from which the lattice artifacts appear to be 0(a 2 ) down to jS = 5.85 and of relatively 
small size. 
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Figure 3: a. Scaling behaviour of m oct ro for 7 fixed values of mpsro (only m^ CAC used), b. Continuum limit 
of /ps^o as function of (mpsro) 2 (only m ? c CAC used). The empty squares are results taken from 1 12| 
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j8 


5.70 


5.85 


6.00 


6.10 


6.20 


6.45 


a (fm) 
L/a 
T/a 


0.171 
12 
32 


0.123 
16 

32 


0.093 
16 

32 


0.079 
20 
40 


0.068 
24 
48 


0.048 
32 
64 


Nconi (me 10 ") 


600 


378 


387 


300 


260 


182 


^conf (m P c CAC ) 


600 


500 


400 




300 





Table 1: Parameters of the simulations 



5. Continuum limit 



Results for the continuum limit of fp$rQ, myr^, m oct ro and nid ec ro are presented in Fig. |3|b 
and |].a. Our determinations of fp$ and my are in good agreement with those from the ALPHA 



Coll. Q12[ ] with non-perturbatively O(a) improved Wilson fermions (the comparison for /ps is 
shown in Fig. |||b). For the chiral extrapolation we used the form fps, my ~ A + Bm 2 PS and m oct , 
tn<iec ~ A + Bmps + Cmp s . In order to compare with the results of the CP-PACS Coll. [ 14 ] , we fixed 
the scale ro through the p mass m p , obtaining ro = 0.576 fm 4 . As a prediction we get f K , rrifj, m\ 
and (working in the SU (3) symmetric limit) fg and m^*. The results can be found in Fig. ||b and 
Tab. U together with those of Ref . [ 14 ] and turn out to be in very good agreement. Notice however 
that in the present work quantities are 0(a) improved and pseudoscalar masses significantly smaller 
than those in Ref. [ |l4| ] (where standard Wilson fermions were used) have been simulated. 



6. Conclusions and Outlook 



In the present study we have extrapolated to the continuum (taking into account possible FSE) 
meson quantities (fp$ and my) and baryon masses (m oct and md ec ) down to pseudoscalar messon 
masses of 235 MeV and 375 MeV respectively. We have presented a strong evidence that lattice 
artifacts are 0(a 2 ) for both definitions of m c (m^ lon and m^ CAC ) and moreover that the use of m^ CAC 
drastically reduces the chirally enhanced 0(a 2 ) lattice artifacts present at small quark masses when 
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U (MeV) 


f K (MeV) 




exp. 


132 


160 


1.22 


tmQCD 


126(5) 


146(3) 


1.15(5) 


CP-PACS 


120(6) 


139(5) 


1.16(3) 



Table 2: Pseudoscalar meson decay constants from the present work and from Ref. jl4| | (in Ref. [[TjJ], 
tadpole-improved perturbation theory has been used to renormalize the axial current). 



using m^ lon . Our results for the continuum extrapolated quantities are in good agreement with 



those from the ALPHA [ ]12[ ] and CP-PACS [ ]14[ ] Collaborations. This is very encouraging in view 
of future dynamical simulations [Q]. There are however other aspects of tmQCD that it is worth to 
investigate, for instance the problem of isospin breaking (see Ref. [|l^] for an exploratory study in 
the quenched approximation). This problem is practically very important for phenomenological ap- 
plications of tmQCD and also strictly related to the phase structure of tmQCD in the neighborhood 
of the critical point []|, [TJ|] and thus directly relevant for dynamical symulations. 
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